Dynamics of colloidal aggregation in microgravity by critical Casimir forces 
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Using the critical Casimir force, we study the attractive-strength dependence of diffusion-limited 
colloidal aggregation in microgravity. By means of near field scattering we measure both the static 
and dynamic structure factor of the aggregates as the aggregation process evolves. The simultaneous 
measurement of both the static and dynamic structure factor under ideal microgravity conditions 
allows us to uniquely determine the ratio of the hydrodynamic and gyration radius as a function 
of the fractal dimension of the aggregate, enabling us to elucidate the internal structure of the 
aggregates as a function of the interaction potential. We find that the mass is evenly distributed in 
all objects with fractal dimension ranging from 2.55 for a shallow to 1.75 for the deepest potential. 
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Since long colloidal aggregation processes have been 
recognized as a class of phenomena that can be described 
by very general rules, often referred to as the universal- 
ity of colloidal aggregation. After the fundamentals in- 
troduced by Verwey and Overbeek [I], the experimental 
discovery that colloidal clusters are endowed with fractal 
structures [2] triggered a lot of work theoretically [3] , by 
simulations 4] and experimentally [5] to describe aggre- 
gation phenomena that are ubiquitous in nature, foods, 
and many other consumer products. Colloidal aggrega- 
tion is central to the formation of gels in systems with 
short-range attraction [6 . However, most of our under- 
standing comes from limiting cases of infinitely strong 
particle attraction, where particles stick irreversibly. The 
current understanding can be roughly summed up by dis- 
tinguishing two different regimes: diffusion limited aggre- 
gation (DLA) and reaction limited aggregation (RLA), 
depending on the slowest phenomenon limiting the clus- 
ter growth |3]. 

One of the big challenges lies in understanding the 
internal structure of the aggregates. Especially at low 
attractive strength, where particles can detach and rear- 
range, the situation is not clear. In this case, the internal 
structure would give important insight into the aggregate 
growth process and the nature of aggregation for a case 
that is most relevant for natural aggregation phenom- 
ena. As shown by Wiltzius [7] , this internal structure can 
be addressed by determining both the gyration radius, 
R g , and the hydrodynamic radius, R^, and their ratio, 
j3, that is directly related to the density-density correla- 
tion function of the aggregate. However, the challenge in 
the determination of f} is that it requires simultaneous 
measurement of both static and dynamic light scatter- 
ing (SLS and DLS) for large objects, which is a difficult 
task. Typically, for such large clusters, the Brownian mo- 



tion becomes very slow and dynamic measurements are 
overshadowed by large-scale convective motion and sed- 
imentation, making the measurement of hydrodynamic 
radii prohibitively difficult. 

Here, we study the internal structure of colloidal ag- 
gregates formed at low to high attractive strength. We 
exploit data obtained on the International Space Station 
(ISS), where there is no convection or sedimentation, and 
pure diffusive motion of the aggregates is guaranteed. 
We realize an effective attractive potential of control- 
lable strength and range by employing critical Casimir 
forces [U [9j . These attractive forces result from the con- 
finement of critical solvent fluctuations between the par- 
ticle surfaces; their strength is determined by the corre- 
lation length of the solvent and hence by temperature. 
This effect thus offers an effective potential that adjusts 
with temperature on a molecular time scale, allowing us 
to study the aggregation process as a function of the 
attractive potential. The application of a recently de- 
veloped method of near field scattering (NFS) provided 
the unique opportunity to have simultaneously both SLS 
and DLS measurements at many angles [10] to deter- 
mine both R g and Rh- This allowed us to elucidate the 
structure of the aggregates under ideal conditions and 
at finite potential. At low attractive strength, particles 
can be expected to detach and rearrange to form denser 
structures. Indeed, from the measured static form fac- 
tors, the fractal dimension df was observed to depend 
upon the interaction strength |llj . Here, we show how 
the particle attraction influences the internal structure of 
the aggregates. We elucidate the internal structure and 
present experimental results for the dependence of the 
ratio /3 = Rh/Rg upon the fractal dimension df. These 
measurements provide a deeper insight into the influence 
of the potential on the aggregation process. 



The experiment, named COLLOID, operated in the 
ESA Mirogravity Science Glovebox, under the ELIPS 
program. Charge stabilized fluorinated latex particles 
400 nm in diameter, with a density of 1.6 g/mL and re- 
fractive index of n p =1.37 [T2] were suspended in a mix- 
ture of 3-methyl pyridine (3MP) in water/heavy water. 
Weight fractions of X^ m = 0.63 for D2O/H2O and X^mp 
= 0.39 for 3MP have been used, the solvent refractive 
index being 1.40. Before adding colloids, the solvent mix- 
ture was purified by distillation under vacuum. Four dif- 
ferent suspensions were prepared, each one with a colloid 
volume fraction of ~ 10~ 4 , and with different salt con- 
centrations of 0.31, 1.5, 2.7 mmol/L of sodium chloride. 
The corresponding Debye screening lengths are, respec- 
tively, Ljj— 14 nm, Ljj= 6.4 nm, Ljj— 4.8 nm. The 
samples were filled into quartz cells under vacuum that 
were tightly sealed. 

We used a collimated laser beam 8 mm in diameter 
with a wavelength of 930 nm to illuminate the sample cell, 
and a 0.25 NA, 20X microscope objective to project the 
transmitted and scattered light onto a CCD sensor. As 
detailed in |llj . near-field scattering images result from 
the interference between the intense transmitted and the 
(fainter) scattered light; these provide directly the scat- 
tered field-field correlation functions of the objects within 
the range of scattering angles selected by the lens. 

We first determined the aggregation temperature, 
T agg , of each sample using the second moment 777,2 of 
the scattered intensity. The second moment is propor- 
tional to the cross-section and the number of scatterers; 
aggregation is detected when the second moment passes 
a threshold of 10~ 3 . We then followed the entire ag- 
gregation process in time after temperature jumps from 
below T agg to T agg , T agg + 0.1, + 0.2, + 0.3 and + OAK, 
increasing the attractive strength with each jump. For 
each temperature, we monitored the aggregation process 
for one hour by acquiring NFS images with a frame rate 
of 1 s _1 in batches of 100. The sample was then cooled 
down to a temperature far below T agg to split up the ag- 
gregates, followed by stirring for at least 3 hours before 
a new measurement was started. 

A typical example of the evolution of 7712 is shown 
in Fig. [T] Here, we have used the reduced time t r — 
(t — to)/t B , where t s is the time a single particle diffuses 
over its own diameter and to is the start of the aggrega- 
tion process, which we define by linear extrapolation of 
the data to vanishing m%. Starting from to, the second 
moment rises continuously, indicating the growth of ag- 
gregates. In this regime, the m 2 curves of all samples at 
all temperatures can be scaled onto a single master curve, 
reflecting the universal growth of the aggregates. As a 
guide and measure of reliability, we also indicate the kur- 
tosis, k = 7714/7712- The kurtosis varies wildly for t < to 
and reaches a defined value of about 3 after that. We 
hypothesize that at these early times before to, subcrit- 
ical nuclei may form and evaporate on the time scale of 
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FIG. 1. Evolution of the second moment and the kurtosis dur- 
ing the aggregation process at the lowest attractive strength 
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FIG. 2. Aggregate growth: Evolution of R g and Rh with time. 
The radius of gyration was obtained from the Fisher-Burford 
fit of the static power spectra. The hydrodynamic radius was 
obtained using eq. [T] 



observation. After this, we can measure a reliable radius 
and diffusion coefficient of the scattering objects. 

Microgravity conditions allow unique measurement of 
the internal structure of the aggregates: they permit the 
slow Brownian motion of large aggregates to be mea- 
sured, until the aggregates exhibit the pronounced form 
factors of fractal objects. The NFS measurement tech- 
nique then allows us to measure the intermediate scat- 
tering function F(q, r) instantaneously with respect to 
the much slower aggregation process and diffusion time, 
enabling us to determine the evolution of both the hy- 
drodynamic and the gyration radius as the aggregates 
grow. 

We first determine the gyration radius from the static 
form factor S(q). This is done using the Fisher-Burford 



fit S(q,R g ) = (1+ {2/3d f )q 2 R 2 g ) df/2 for fractal aggre- 
gates that depends only on the radius of gyration and the 
fractal dimension, df |llj . The resulting evolution of R g 
is shown by open symbols in Fig. [2J In good approxima- 
tion, R g grows as a power law, R g ~ t r s , as expected 
for a pure DLA process. In pure DLA, the inverse expo- 
nent equals the fractal dimension of the aggregate. We 
find that, indeed, the slopes in Fig. [2] are in good agree- 
ment with the fractal dimension df — 2.4 and df = 1.8 
determined from the static structure factor respectively 
at T agg and T agg + OAK [TT]. Hence, we can describe 
the entire aggregation processes within the framework of 
DLA with one consistent value of the fractal dimension. 
To determine the hydrodynamic radius, we relate the 
decay of the intermediate scattering function F(q, t) to 
the diffusion coefficient, D. In the case of monodisperse 
simple spherical objects, D is related to the hydrody- 
namic radius, Rh, through the Stokes-Einstein relation 



D = 



k b T 



in which r\ is the solvent viscosity. No depen- 
dence of Rh on q is expected. For our aggregates, how- 
ever, Rh as obtained from the decay of F(q, t) varies with 
q: the additional rotational degrees of freedom of objects 
with internal structure lead to a faster decay of F(q,t). 
In this case, the effective diffusion coefficient -D e // a s 
defined from the decay time r via r = \/D e ffq 2 is no 
longer related to the hydrodynamic radius via the Stokes- 
Einstein relation. To account for the rotations quanti- 
tatively, we define a class of spherical symmetries that 
determines how much an object needs to rotate before 
it decorrelates. The higher the symmetry, the smaller is 
the angle and thus the time needed to decorrelate the 
field [13]. Because the symmetry is connected to qR g 
via the scale-invariant fractal structure, qR g determines 
uniquely the correction to the decorrelation times. This 
allows us to relate -D e // to D via the static structure fac- 
tor S(q) [H]. For S(q), we use the structure factor as 
measured at the same time. To be independent of fluc- 
tuations in S(q), we again use the Fisher-Burford fit as 
detailed above; this expression provides a good fit to the 
measured average structure factor of the aggregates. Us- 
ing it, we can rewrite eq. 7 from Lin et al. [14] into the 
following form: 
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where we have separated static and dynamic quanti- 
ties. The right-hand side is given by the static data; it 
provides a master curve solely dependent on the product 
qRg, taking into account both translations and rotations 
of the aggregates. The left-hand side is determined from 
the dynamics: taking R g from the static data, the only 
free parameter is the hydrodynamic radius. To determine 
it, we define f(Rh) that equals the right and left-hand 
side of the equation, and determine diffusion coefficients 
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FIG. 3. Master curve according to eq. [T] of the scattering 
of aggregates growing by diffusion-limited aggregation at the 



lowest attractive strength at T = T a 



D e ff from the time decay of F for 60 different g-values 
in the range 0.8 - 1.9 /jm" 1 . We then adjust Rh so that 
f{Rh) from the dynamic data provides the best fit to the 
master curve. 

An example of the dynamic data compared to the mas- 
ter curve is shown in Fig. [3J The weakest interaction is 
considered, namely at the temperature T = T agg . The 
data is well fitted with the ratio Rh/R g = 1.05, one con- 
sistent value for the entire aggregation process. Fits of 
similar quality are obtained for all other interaction po- 
tentials as long as the aggregates are not too large. The 
situation changes for the highest interaction potential, 
when aggregates grow to larger sizes: the deviation from 
the master curve becomes larger due to polydispersity 
effects [7 J. By taking the polydispersity into account we 
can reduce the deviations and again obtain a good fit for 
one value of /?. We accounted for polydispersity explicitly 
assuming a cluster mass distribution of the form: 



N(M) = 



N T 

W) 
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where N(M) is the number of clusters of mass M, Nt is 
the total number of clusters, and (M) is the average clus- 
ter mass, which is related to the radius of gyration R g via 
(M) = (R g /a) df , where a is the monomer radius. This 
distribution represents a good description for the DLA 
case (see [2] and references therein for details). To de- 
termine the effective diffusion coefficient, for each cluster 
mass, we determined the individual diffusion coefficients 
using eq. [TJ and we averaged the diffusion coefficients 
of all clusters, taking into account their statistical weight 
given by N(M) [T3J[T5]- The results of this method agree 
with the earlier one that neglects polydispersity to within 




FIG. 4. Ratio of hydrodynamic to gyration radius as a func- 
tion of df for three different salt concentrations (see text for 
details). Lines indicate the expected dependence obtained for 
unit step, Gaussian and exponential density-density correla- 
tion functions (from top to bottom). Insets show holographic 
reconstructions of aggregates grown at T — T agg +0AK (high- 
est attraction, top) and T = T agg (lowest attraction, bottom). 
The length of the scale bar is 25^m. 



2-3% for small values of qR g ; however, by incorporating 
polydispersity, we now obtain a similarly accurate fit over 
the entire range of qR g , even for large values of qR g . We 
note that this large-aggregate regime is important in mi- 
crogravity measurements such as the one explored here, 
because large aggregates do not settle on the observation 
time scale and stay in the field of view. In contrast, this 
regime is rarely approached on earth, and is not reached 
in [H] because of the fast settling of large aggregates. 

Thus, we have obtained the hydrodynamic radius upon 
growth of the aggregates for a large range of sizes. The 
resulting values of Rh as a function of time are indicated 
by closed symbols in Fig. [2] Comparison with R g in the 
same figure shows that the ratio Rh/R g becomes con- 
stant for larger aggregates, being /3 = 1.1 ± 0.012 for 
T = T agg , and /3 = 0.85 ±0.015 for T = T agg +0AK. For 
smaller aggregates, the ratio Rh/R g appears not constant 
because the aggregates cannot be considered spherical as 
the monomers are 0.4 \im in diameter and therefore ag- 
gregates consist of a few monomers only. This small size 
effect is also observed in the growth as a deviation from 
the expected power law R g = at 1 ^ 1 dependence. By 
considering only sufficiently large aggregates that have 
a well-defined fractal dimension and structure, we can 
now elucidate the internal structure as a function of the 
attractive potential and fractal dimension. 



We plot the ratio 



Rh/R g as a function of d 



in Fig. HI Measurements at all temperatures from T a 



I 



2.7mmol/L, green triangles) are included. A systematic 
dependence of j3 on df is observed, indicating a consis- 
tent internal structure independent of the salt concen- 
tration. To interpret the /3-values, we indicate by lines 
the dependence for a fully compact object with a unit- 
step density distribution (solid curve), for a fluffier object 
with a Gaussian density profile (dashed curve) , and for an 
exponential profile (dotted curve) [7j . The data lies clos- 
est to the unit step, indicating that the aggregates have 
fairly compact internal structure, regardless of the fractal 
dimension and thus the attractive strength. To illustrate 
the shape of the aggregates, as insets we show results 
obtained by holographic reconstruction of the aggregates 
at the latest stages of aggregation for the two conditions 
considered in Fig. [2] The fluffier structure of aggregates 
formed at higher attractive interaction strength is clearly 
observable. 

In summary, we measured the internal structure of 
DLA clusters over a wide range of attractive potential 
strength. This was possible owing to the peculiarity of 
the critical Casimir effect allowing us to induce tunable 
interactions that result in aggregates with a wide range 
of fractal dimensions. Notice that no systematic depen- 
dence on the added salt is observed, and the only depen- 
dence is on the fractal dimension df, set by the depth of 
the potential well [TT]. The salt concentration i.e. the 
repulsive part of the potential has no effect on df , as it is 
only the depth of the potential well that determines the 
fractal dimension. While the aggregation temperature is 
set by the ratio of the solvent correlation length to the 
Debye screening length, it is the temperature increment 
from T agg that sets the depth of the potential well, and 
consequently determines the aggregate structure, i.e. its 
fractal dimension. This is the first time that a DLA pro- 
cess is studied with the widest range of possible fractal 
structures, owing to the variable critical Casimir poten- 
tial, microgravity and the near field scattering technique 
that measures the intermediate structure function as a 
function of time. 
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